read_data;

%% set parameters

% externally calibrated parameters
R=1/1.02;
beta=1/R;
T=12; % duration of investment period

% estimated technology parameters (Table 5, kappa=1)
rho=1-1/0.2;
gma=1-1/0.5;

delta1=0.17;
delta2=0.83;

phi_m=[8.2 0.1 -0.96 -2.42 -0.5 -1.7 -0.59 0.46]';
phi_f=[4.11 -0.47 0.64 -1.04 -0.95]';
phi_Y=[-1.19 0.61 0.08 0.1 0.01 -0.23 -0.07 0.06 0 -0.63]';
phi_theta=[1.15-0.3 -0.2 -0.04 -0.42 0.05 -0.08 -0.14 0.06 -0.01 0.17]';

% preference parameters to calibrate
alpha=zeros(N,1);
psi_m=zeros(N,1);
psi_f=zeros(N,1);
share_model_age=zeros(2,2,T,3);

%% start running
set_regressors;

if calibration
    calibrate;
else
    read_parameters;
end    

if decomposition
    decompose;
elseif price_reduction
    reduce_price;
elseif free_childcare
    solve_free_childcare;
end